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Abstract 

We estimate character sums with n!, on average, and individuahy. 
These bounds are used to derive new resuhs about various congruences 
modulo a prime p and obtain new information about the spacings be- 
tween quadratic nonresidues modulo p. In particular, we show that 
there exists a positive integer n <C p^/'^^^, such that n! is a prim- 
itive root modulo p. We also show that every nonzero congruence 
class a ^ (mod p) can be represented as a product of 7 factorials, 
a = ni\...n-j\ (modp), where maxjnj | i = 1,...7} = 0(j?-^^/^^+^), 
and we find the asymptotic formula for the number of such represen- 
tations. Finally, we show that products of 4 factorials ni\n2\n^\rn\, 
with maxjni, 712, ^3, = 0{p^/'^^^) represent "almost all" residue 
classes modulo p, and that products of 3 factorials ni\n2\n^\ with 
max{ni, 712, 71-3 } = 0(p^/^+^) are uniformly distributed modulo p. 
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1 Introduction 

Throughout this paper, p is an odd prime. Very little seems to be known 
about the distribution of n! modulo p. In Fll in 0, it is conjectured that 
about p/e of the residue classes a (mod p) are missed by the sequence n\. 
If this were so, the sequence n! modulo p should assume about (1 — l/e)p 
distinct values. Some results in this spirit appear in 2J. The above conjecture 
immediately implies that every residue class a modulo p can be represented 
as a product of at most two factorials. Unconditionally, it is easy to see that 
three factorials suffice. Indeed, = p\ (mod p), and, as it has been remarked 
in [S], equation (5), see also (2j, the Wilson theorem implies that 

b\ ■ [p - 1 - b)\ = {-if^^ (modp) (1) 

holds for any b E {1, ... ,p — 1}. Therefore, if a G {1, . . . ,J9 — 1}, then with 
b = a^^ (mod p), we have 

a={{p-iy.Y'{b-l)\-{p-l-by. (modp), (2) 

where r;, G {0,1} is such that = b + 1 (mod 2). However, the above 
argument does not apply to proving the existence of representations involving 
factorials of integers of restricted size, neither can it be used for estimation 
of the number of representations. 

In this paper, we first estimate character sums with n\ on the average, 
and individually. We use these estimates to show that for every e and p 
sufficiently large, there exists a value of n with n = 0{p^^'^^^) and such that 
n\ is a primitive root modulo p. 

We apply these estimates to prove that every residue class a ^ (mod p), 
can be represented as a product of 7 factorials, a = nil . . .nj\ (mod p) with 
max{nj | i = 1, . . . , 7} -C pii/^^+e^ jf only want that "most" of the residue 
classes modulo p be represented as a product of factorials in the same range 
as above (and even a slightly better one), then we show that four factorials 
suffice. Moreover, our results imply that for every e > and sufficiently 
large p, every residue class a ^ (mod p) can be represented as a product 
of i = \_e~^\ + 5 factorials, a = nil... nil (modp), where maxjnj | i = 
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1, ...,£} <^ g^jgQ show that products of three factorials ni!n2!n3!, 

with max{ni, n2, Tis} = 0{p^^^~^'^), are uniformly distributed modulo p. 

Our basic tools are the Weil bound for character sums, see [T^IT^ ITH]. and 
the Lagrange theorem bounding the number of zeros of a non-zero polynomial 
over a field. 

Some of the results of this paper have found applications to the study 
of arithmetic properties of expressions of the form n! + f{n), where f{n) 
is a polynomial with integer coefficients (see fl]), or a linearly recurrent 
sequence of integers (see ^H])- In particular, an improvement of a result of 
Erdos and Stewart [S], obtained in [T3|, is based on these results. 

Throughout the paper the implied constants in symbols 'O' and '-C' may 
occasionally, where obvious, depend on integer parameters i and d and a 
small real parameter e > 0, and are absolute otherwise (we recall that A <^ B 
is equivalent to A = 0(B)). 

Acknowledgements. The authors would like to thank Vsevolod Lev 
for several useful comments. During the preparation of this paper, F. L. was 
supported in part by grants SEP-CONACYT 37259-E and 37260-E, and I. S. 
was supported in part by ARC grant DP0211459. 



2 Character Sums 

Let Wp be a finite field of p elements. We always assume that Wp is repre- 
sented by the elements of the set {0, 1, ... ,p — 1}. 

Let X denote the set of multiplicative characters of the multiplicative 
group F* and let X* = X\{xo} be the set of nonprincipal characters. 

We also define 

e{z) = exp{2TTiz/p), 

which is an additive character of Wp. 
It is useful to recall the identities 



^X{u) 



and 

p-i 



e[au) 

a=0 



0, if u ^ 1 (mod p), 

p — 1, if u = 1 (mod p), 



0, if M ^ (mod p), 
p, if M = (mod p), 



3 



which we will repeatedly use, in particular to relate the number of solutions 

of various congruences and character sums. 

Given X G A", a polynomial / G Fp[X], and an element a G Fp, we 
consider character sums 

H+N 

T{xJ,H,N)= J2 X(n!)e(/(n)) 

n=H+l 

where we simply write T{x, H, N) is / is identical to zero, and 

H+N 

S{a, H,N)^ e (an!) . 

n=H+l 

We obtain a nontrivial upper bound for "individual" sums T{x, f, H, N), 
and also nontrivial upper bounds for the moments of T(x, /, H, N) and 
S{a,H, N). 

Theorem 1. Let H and N be integers with 0<H<H + N<p. Then for 
any fixed integer d> 0, the following bound holds: 

max max \T{x, /, H, N)\ « N'/y/%logpy/\ 

deg f=d 

Proof. For any integer /c > we have 

H+N 

T{x,f,H,N)^ J2 x{{n + k)\)e{f{n + k)) + Oik). 

n=H+l 

Therefore, for any integer K > 0, 

T{x,f,H,N) = ^W + 0{K), (3) 



where 



K-l H+N 



^ = E E X{{n + k)\)e{f{n + k)) 

k=0 n=H+l 

H+N K / k \ 

= E E^h'n(^+^)p(/(^+^)) 

n=H+l k=l \ i=l J 
H+N K-l / k \ 



n=H+l k=0 \i=l 
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Applying the Cauchy inequality, we derive 

H+N 

\W\'^ < N 

n=H+l 
K-1 H+N 

= 5^ \(*fc,™H)e(/(n + A;)-/(n + m)), 

k,m=On=H+l 

where 

k m 

= l[{x + i)l[{x + j)-' 
i=i j=i 

and hereafter * means that the poles of the corresponding rational function 
are excluded from the summation (we also recall that |zp = zz for any 
complex number z, and that x{(^) = x{(^^^) holds for every integer a ^ 
(mod p) where x is the conjugate character of x)- 

Clearly, ii K < p then, unless k = m, the rational function \&fc ,„(X), has 
at least one single root or pole, and thus is not a power of any other rational 
function modulo p. 

For the 0{K) choices ofO<A; = m<i^' — 1, we estimate the sum over 
n trivially as N. 

For the other 0{K'^) choices ofO<k,m<K — 1, using the Weil bound, 
given in Example 12 of Appendix 5 of |18] (see also Theorem 3 of Chapter 6 
in or Theorem 5.41 and the comments to Chapter 5 of Jl]), we see that, 
because x ^ then, for any a e Fp, we have 

Y,*^ (*fc,m(^)) e(/(n + k)- f{n + m) + an) = 0{Kp^/^). 

n=0 

Therefore, by the standard reduction of incomplete sums to complete ones, 
(see P), we deduce 

H+N 

Y *^ (^fc,m(n)) e(/(n + k) - f{n + m)) = 0{Kp^''^ \ogp). 

n=H+l 

Putting everything together, we get 

<t:N [KN + K^p^'"^ \ogp) . 



K-l 



k=0 



,i=l 
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Therefore, by 0, we derive 



Taking K = [A^^/^p ^^'^{logp) ""^^^J, we finish the proof. 



□ 



It is clear that for any e > there exists some 5 > 0, such that if 



provided that p is large enough. 

Clearly, Theorem Q immediately implies that among the values of n\, 
where n = H + 1, . . . , H+N, there are N/2 + [N^/^p^/^{\ogpy/'^) quadratic 
residues and nonresidues. Remarking that each change in the value of the 
Legendre symbol {n\/p) corresponds to a quadratic non- residue n we can 
derive a certain result about the distribution of spacings between quadratic 
non-residues modulo n which does not seem to follow from any of the previ- 
ously known results, see [TT] . 

Let rij be the jth quadratic nonresidue modulo p and let dj = rij — rij^i, 
the jth spacing, j = 1, . . . , (p — l)/2, where we put uq = 0. 

Corollary 2. Let J be an integer with p^^'^logp < J < {p — l)/2. Then the 
following bound holds: 



N > then 



\T{xJ,H,N)\<Np 



-s 



J 



Y,{~^yd,<^.f'V'\\ogpf'\ 



Proof. We have 




From the Polya 



Vinogradov hound 





we see that nj 
result. 



2 J + 0{j> 



log p) <^ J and by Theorem ^ we derive the 

□ 
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Obviously 

J 

Y^d,=nj = 2J + 0{p^/Hogp), 
i=i 

which demonstrates that for every J > p^l'^+^ the odd and even spacings dj, 
j = 1, . . . , J, are of approximately the same total length. 

We now denote by Q{H, N) the number ofn = H + l,...,H + N such 
that n! is a primitive root modulo p. 

Corollary 3. Let H and N be integers with 0<H<H + N<p. Then, for 
any fixed e > 0, the following bound holds: 

Q{H, N) = N^^^P^ + O (iV3/V/8+-) . 
p — 1 

Theorem 4. Let H and N be integers with 0<H<H + N<p. Then for 
any fixed integer d > 0, the following bound holds: 

max V |r(x, /, H, N)\' « pN'/'. 

deg/=d ^ 

Proof. Arguing as in the proof of Theorem ^ and applying the Holder in- 
equality to 0, we derive that for any K 

K-l H+N 

\T{x,f,H,N)\' « K~'nJ2 

k,m=0 n=H+l 

e{f{n + k)- f{n + m)) + K\ 

where 

k m 

1=1 j=l 

and as before ^ * means that the poles of the corresponding rational function 
are excluded from the summation. Therefore, 

\T{x,f,H,N)\' 

K-l H+N 

« K-'N J2 E * E ^ i^kAn)) + PK\ 

k,m=On=H+l xeA' 
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The sum over x vanishes, unless 

^K^in) = 1 (mod p) (5) 

in which case it is equal to p — 1. For K pair {k, m) with k = m then there 
are N possible solutions to for other O^K"^) pairs there are at most K 
solutions to (0). Thus 

J2\nx,f,H,N)\' « K'^N{K^ + KN)p + pK' 

= {NK + N^R-^ + K^) p. 
Taking K = [A^^/^J , we finish the proof. □ 

Theorem 5. Let H and N he integers with 0<H<H + N<p. Then for 
any fixed integer i > 1, the following bound holds: 

p-i 

J2\S{a,H,N)\'<^pN'/'. 

a=0 

Proof. Arguing as in the proof of Theorem ^ we derive 

p-l K-l H+N p-l 

J2\S{a,H,N)\'<t:K-'N E E ^ ("^'^'^-'-H) + P^'' 

a=0 k,m=0 n=H+l a=0 

where 

i=i j=i 

The sum over a vanishes, unless 

n\^k,m{n) = (modp), (6) 

in which case it is equal to p. 

As before, we see that is a nonconstant polynomial of degree 

0{K), unless k = m Because n! ^ (mod p) for 0<H<n<H + N<p, 
we derive 

p-i 

Y,\S{a,H,N)\^ <t: {NK + N^K-' + K^)p. 

Taking K = [A^^/^J and remarking that with this value of K the last term 
never dominates, we finish the proof. □ 
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3 Sums and Products of Factorials 



For integer i > 1 and H and with < H < H + N < p let ns denote by 
l£{H, N) and Jii^H, N) the number of solutions to the congruences 



21 



Y[ni\=Y[ni\ (modp), H + 1 < rii, . . . ,n2e < H + N, 
1=1 i=i+i 



and 

i 2i 



rij! = ^ rij! (mod p), if + 1 < rii, . . . , ?t,2£ < if + iV, 



i=\ i=e+i 



respectively. 

From the properties of multiplicative and additive characters we imme- 
diately conclude that 

and 



-J2\S{a,H,N)f' = MH,N). (8) 

^ a=0 

The same arguments as in the proof of Theorem El lead to the bound 

We now show that for Ii{H, N) one can derive a more precise estimate. 

Theorem 6. Let H and N be integers with 0<H<H + N<p. Then for 
any fixed integer i > 1, the following bound holds: 



Proof. We prove this bound by induction. If £ = 1 then Theorem E] taken 
with f{X) = 0, together with ((Zj) immediately imply the desired bound 
Ii{H,N) < A^3/2_ 
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2^-3+2- 



Now assume that £ > 2 and that h-i{H, N) < A^' 
K < N and note that by the Cauchy inequahty we have 



. We fix some 



H+N 


2 






n=H+l 









K 



k=l H+(k-l)N / K<m<H+kN / K 



K 



k=l 



H+{k-l)N/K<m<H+kN/K 



Therefore 



Ie{H,N) 



K 



P 



H+{k-l)N/K<m<H+kN/K 



H+N 
n=H+l 



X{n\ 



k=i xeA- 
= KIe{K,H,N), 

where Ie{K, H, N) is the number of solutions to the congruence 

£-1 2£-2 

mi\Y\_ni\ = Y\_ ^i^- (mod p) 



2^-2 



1=1 



with + 1 < ni, . . . , n2£-2 < H + N and H + {k ~ 1)N/K < mi,m2 < H + 
kN/K for some k = 1, . . . , K. For each of N pairs (mi, with mi = m2, 
there are exactly /^_i(if, A^) solutions. Also we see that if ni, . . . , ^2^-2 are 
given then for each fixed value of r = mi — m2, there are no more then 
\r\ values solutions in mi,m2 (because at least one of mi, 1712 satisfies a 
nontrivial polynomial congruence of degree |r|). Certainly r = 0{N/H). 
Putting everything together and using the induction assumption we obtain 



h{K, H, N) < NI,^i{H, N) + {N/KfN' 



21-2 



N' 



2^-2+2-^+1 _^ jY2^_^-2^ 



Therefore Ie{H, N) < kN^^-^+^' 
we obtain the desired bound. 



N^^K-\ Choosing K 



N 



l-2~ 



n 



We now show that, for A^ > p^/2+^ the above bound on l£{H,N), com- 
bined with Theorem produces an asymptotic formula for Ii{H,N). In 
particular for /7 = 0, A^ = p— 1, this asymptotic formula is nontrivial for 
£ > 4. 



10 



Theorem 7. Let H and N be integers with 0<H<H + N<p. Then for 
any fixed integers i > r > 1, the following bound holds: 



I^^H, N) = ^ + O { N^'l^+'^l^'^ 
p —1 V 



Proof. Similar to 8j, we have 



+2-'-p(^-r)/4(i^gp)(£-r)/2 



V — I ^-^ 



p — i p 



H+N 



n=H+l 



21 



- 1 ^ n- 1 ^ 



x&x* 



H+N 
n=H+l 



21 



N' 



21 



P — 1 



+ 



( 1 


H+N 


2£-2r 


H+N 


2r\ 


max 

I p - 1 X6^* 


n=H+l 


E 


n=H+l 





(note that in the last sum we bring back the term corresponding to x = Xo)- 
The result follows from ((Zj) and Theorems C] and IHl □ 

In particular, using Theorem [7| with r = 1 we obtain 

h{0,p-l)=p' {l + Oip-'/\\ogpf')) . 
We now denote by Fi{a, H, N) the number of solutions to the congruence 



JJ^nj! = a (modp), H + 1 < ni, . . . ,ng < H + N, 



i=l 



where a G F*. 



The same arguments as the ones used in the proof of Theorem [7| imply: 

Theorem 8. Let H and N be integers with 0<H<H + N<p. Then, for 
any fixed integers i >2r > 1, the following bound holds: 



max 

a&W* 



F,{a,H,N)- 



p — 1 



^ JY3£/4+r/2-l+2-'-^(^-2r-)/8QQgp^)(^-2r)/4_ 
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In particular, using Theorem |H1 with r = 1 we obtain 

F,ia, 0,p - 1) = /(I + O {p-'/\\ogpf') , 
and for any e > 0, using Theorem |H1 with r = 2 we obtain 

Fr{a, H,N) = -—{1 + o(l)), for > 
P 

hold for all a G F*. 

Let Vi{H, N) be the number of a G F* for which Fe{a, H, N) > 0, that 

is, 

T4(i/,iV) = #|j]nJ (modp), \ H + 1 < n,, . . . ,n, < H + 

Theorem 9. Let H and N be integers with 0<H<H + N<p. Then for 
any fixed integers i > r > 1, the following bound holds: 

Ve{H, N)=p + (^iV-^/2+r/2-i+2-'-^(^-.+8)/4^i^gp)(^-0/2j 
Proof. We may assume that i > 2, otherwise there is nothing to prove. Let 

£ = |/i G Fp I ^ JJnJ (modp), H + 1 < n^, . . . ,ne < H + . 

Then, 

^ xeX ni,...,ne=H+l he£ 

Separating the term corresponding to Xo and, for x ^ applying Theo- 
rem [T] to the sums over ui, . . . , ni^r, we obtain 



4 = 1 



H+N 



p — 1 



< (iv^/y/siogV^p 



t-r \ 



p 



H+N 



E 



n=H+l 



As before, we now extend summation over all characters x G A" and by the 
Cauchy inequality, we derive from (jZj) and Theorem IHl 



E^(^) 



H+N 

E 

n=H+l 



r\ 2 



E^(^"; 



hee 



H+N 



E 

n=H+l 



2r 



2 Ar2r-l+2- 
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Therefore, 



'v'#^.<fiv»/V'»iog'/v^'"' 



P- 

which finishes the proof. 



2r-l+2- 



1/2 



□ 



In particular, using Theorem El with r = 2 we see that for > p6/7+£^ 
we have that only o(p) residue classes modulo p cannot be represented as 
ni\n2\n^\ni^\ (mod p) with H + 1 < ni, n2, ns, < H + N. 

We recall that the discrepancy D of a sequence of M points {'yj)fLi of the 
unit interval [0, 1] is defined as 



D 



sup 



A{1) 



M 



1X1 



where the supremum is taken over the interval I = [a,f3] C [0, 1] of length 
\I\ = P — a and A{X) is the number of points of this set which belong to X 

(see mini). 

For an integer a with gcd{a,p) = 1, we denote by De{a, H, N) the dis- 
crepancy of the sequence of fractional parts 



P 



n 

i=l 



H + 1 <ni,...,ni< H + N. 



Obviously, 

D,{a,H,N) 



max 

0<h<k<p-l 



k 

J2 Fe{a~'c,H,N)- 



k-h 



c=h+l 



P 



+ Oip-'), (9) 



thus Theorem IHl can be used to estimate Di{a, H, N). However, we show that 
the Polya-Vinogradov bound (jH) leads to stronger results. 

Theorem 10. Let H and N he integers with 0<H<H + N<p. Then for 
any fixed integers i >2r > 1, the following bound holds: 



max \De{a, H,N)\<^ p(i-2r+4)/8^^^^pY'- 

l<a<p— 1 



2r+4)/4 
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Proof. We have 



J2 Fe{a-'c,H,N)- 



k-h 



c=h+l 



P 



^ k H+N I 

(Ay EE E x(«.-n".! 



^ k f H+N 

E ^(«) E E 

XeA"* a=h+l \n=H+l 



i=l 



H+N 


e-2r 


H+N 


2r 




E 






n=H+l 




n=H+l 





Thus, applying the bound (jH), we deduce 
J2 F,{a-h,H,N)-'^N^ 

c=h+l ^ 



<^ max 



and the result follows from ((Tj) and Theorems ^ and IHl 

In particular, using Theorem with r = 1 we obtain that 
max |D3(a,0,p- 1)1 = O {p~^'^{\ogpf'^) , 

l<a<p— 1 

and also that for any £ > 0, 

max |L)3(0,p- 1)1 = o(l), for > ^5/^+^ 

l<a<p— 1 



We also note that Theorem ITOl implies that 



max 

l<a<p-l 



H+N 



n=H+l \ i=l 

for £ > 2r > 1. 



□ 



2r+4)/4 



Let G£(a, A^) be the number of solutions to the congruence: 

JJ^ = a (mod p). 
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in positive integers ni, . . . , with 



i=l 



It has been shown in [Tn| that for any e and sufficiently large p, Ge{a, N) > 
provided that i > and N - i > p^/^+\ In 0, the same result has been 
obtained under a much weaker condition N—i > p^/^+^ , Here, concentrate on 
the value of £ and show that it can be taken as £ = 0(1) provided > p^/^+s. 

Theorem 11. For any fixed integer £ > 1 and any integer N with 1 < N < 
p/i, the following bound holds: 



max 
aeip; 



G,{a,N) 



N -1 



p-l\i-l 
Proof. For a ^ (mod p), we have 



< A^3^/4p(^+6)/8QQgp)(^-2)/4_ 



Giia,N) 



p 



ni,...,nf>l y£X \ i=l / 



ni,...,nt>l X^X 
n\ + ...-\-ni='N 



where the sum is taken over all multiplicative characters x modulo p. Sepa- 
rating the contribution from the principal character Xqi we obtain 



G,(a, iV)- 



1 (N -\ 



p-\\i-\ 



< 



1 



p — 1 



where 



R 



ni,...,ni>l \j=l 
ni+...+nt=s 



p-1 



XeX* ni,...,ni=l \i=l J ^ c=0 

because if iN < p then the congruence ni + . . . + ue = s (mod p) with 
1 < ni, . . . ,n£ < N is equivalent to the equation ni + . . . + n£ = s. Therefore, 

e 



1 



c=0 x^X* 



^Xinl) e{cn) 



n=l 
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Arguing as in the proof of Theorem [3 we derive the result follows from ((Tj) 
and Theorems Q and IHl □ 

For example, for any fixed e > and p/i > N > have 

G.(a,iV) = ^('Jj^') (1 + 0(1)) 

for every fixed £ > + 4. 

We remark that one can easily drop the condition N < p/£ in Theoremllll 
Let F{a, H, N) = Fi{a, H, N) be the number of solutions of the congru- 
ence n\ = a (mod p),H+l<n<H + N. 

Theorem 12. Let H and N be integers with 0<H<H + N<p. Then 
following bound holds: 

max F (a, H,N) < iV^/^. 
aeip; 

Proof. Let > be a parameter to be chosen later. Let 

A={H+l<n<H + N\ a = n\ (mod p)} = AiU A2, 

where 

Ai = {n E A \ \n ~ m\ > K for all m 7^ m G A} and A2 = ^\.4,i. 

It is clear that i^Ai <^ N/K. Assume now that n G ^2- Then there exists 
a nonzero integer k with \k\ < K and such that n\ = {n + k)\ (mod p). For 
each fc, the above relation leads to a polynomial congruence in n of degree 
\k\ and therefore it has at most \k\ < K solutions n. Summing up over all 
values of k with \k\ < get that #v42 -C K'^. Thus, 

F(a, i/, AT) = #^ = #^1 + < ^ + K\ 

and choosing = [A^^^'^J we get the desired inequality. □ 

We have seen that the Wilson theorem immediately implies the inequality 
^2(0, p — 1) > {p ~ l)/2- We now show that this bound can be slightly 
improved. 
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Theorem 13. The following bound holds: 



V2{0,p-1) > -p + 0{p 



Proof. By we see that if a = 6 ^ (mod p),l<b<p— lis odd then 
a G V2{0,p-l). 



By (Q, we see that if a = c~^(c + 1)^^ (mod p) with some even c = 2u, 
1 < c < p — 3, then a G V2(0, p — 1) too. Thus each such c which corresponds 
to an even 6 = 2f in the above representation, contributes one new element 
to V2{0,p — 1). It is also clear that no more than two distinct values of c can 
contribute the same element. 

Therefore, ¥2(0, p — 1) > (p — l)/2 + W/2, where W is the number of 
solutions of the congruence 



The Weil bound yields W = p/A + 0{p^^'^\o^ p) (see |1J), which concludes 



We remark that Theorem immediately implies that for every integer a 
there exists a representation a = ni\n2\ + n^lnj (mod p) with some positive 
integers ni, n2, n^, n^. 

4 Concluding Remarks 

Most of our results hold in more general settings. For example, let m > 1 be 
any fixed positive integer and put 



Then Theorem ^ holds with T{x, f, H, N) replaced by T(m, x, f, H, N). In 
particular, if we write Q{m, H, N) for the number ofn = H + l,...,H + N 
such that n\, . . . , [n + m — 1)\ are all primitive roots modulo p, then the 
estimate 



2u{2u + 1) = 2v (mod p) 



0<u,v< (p-3)/2. 



the proof. 



□ 
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holds for any fixed integer i > 1. 

Let Q be the set of all distinct prime divisors of p — 1. For a set 7?. C Q, 
we denote by T{TZ, H, N) the number oin = H + l,...,H + N such that for 
every g G Q, n! is a gth power residue modulo p if and only if g G 7?.. Then 
the estimate 

T{n,H,N) = Nl[^ H i + o(iVi-i/V+2)/4£(m)) 
qen ^ qeQ\n ^ 

holds for any fixed integer i > 1. 

Techniques of the present paper apply also to the sequences 

2n\ (2n)\ , , , 

^ - ^ ' {2n+l)\\ = 1 - 3... ■ (2n + 1), 



n J {n\y ' 

and many others, as well as their combinations. 

Also, with some minor adjustments, our methods can be used to obtain 
similar, albeit somewhat weaker results for composite moduli. In this setup, 
our basic tools such as the Weil bound and the Lagrange theorem, have to be 
replaced with their analogues in residue rings modulo a composite number. 
See, for example, P| for bounds of character sums, and [HI for bounds on the 
number of small solutions of polynomial congruences. 

While the results of the present paper represent some progress towards 
better understanding the behaviour of n\ modulo p, there are several chal- 
lenging questions that deserve further investigation. For example, our Theo- 
rem IT2l gives a nontrivial upper bound on F{a, H, N), but we conjecture that 
this result is far from being sharp. We do not have any nontrivial individual 
upper bounds for S{a, H, N). 

Certainly, studying Vi{H,N) is of primal interest. Trivially, we have 
Vi{H, N) > (A^ — 1)^/^ (to see this it is enough to recall that n = n\/{n~l)\), 
but we have not been able to obtain any better lower bound. In the opposite 
direction, answering a question of Erdos, Rokowska and Schinzel jT7] have 
showed that if the residues of 2!, 3 !,..., (p — 1)! modulo p are all distinct, 
then the missing residue must be that of — ((p — l)/2)!, that p = 5 (mod 8), 
and that no such p exists in the interval [7, 1000], but it does not seem to 
be even known that there can be only finitely many such p, or, equivalently, 
that Vi(0,p — 1) = p — 2 can happen only for finitely many values of the 
prime p. 
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It is very tempting to try to generalize the proof of Theorem and 
consider longer products c(c+l) . . . (c+m). This may lead to an improvement 
of the constant 5/8 of Theorem IT^ However, to implement this strategy one 
has to study in detail image sets of such polynomials (and their overlaps), 
which may involve rather complicated machinery. 
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